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Abstract 

We describe the meson-meson data for the {IJ^^ = 00++) wave at 280 < ^/s < 1900 
MeV in two approaches: (i) the K-matrix approach and (ii) the dispersion relation D- 
matrix method. With a good description of low energy data (at 280 < ^/s < 900 MeV) as 
well as the data of two-meson transition amplitudes and antiproton-proton annihilation 
into three pseudoscalar meson states (at 450 < ^/s < 1950 MeV) we have found the 
positions of the resonance poles: (i) for the a meson pole: M^- = (390 ± 35) — i(235 ± 50) 
MeV; (ii) two poles for the /o(980), on the second sheet (under the vrvr cut): Mj = 
(1011 lb 5) — i{35 lb 5) MeV, and on the third sheet (under the tttt and KK cuts), Mjj = 
(1035 lb 50) - i(460 lb 50) MeV; for the /o(1370) meson, M = (1285 ib 30) - i(160 ib 20) 
MeV; for the /o(1500) meson, M = (1488 ib 4) - i(53 ib 5) MeV; for the /o(1790) meson, 
M = (1775 lb 25) - i(UO ± 15) MeV and for the broad state /o(1200 - 1600) M = 
(1540 lb 120) — i(550 lb 70) MeV. Our estimation of the scalar-isoscalar scattering length 
obtained under different parameterizations and assumptions about the quality of low 
energy tttt scattering data is Oq = (0.215 =b 0.040)/i^^. We also discuss the idea according 
to which the o"- meson could be a remnant of the confinement singularity, 1/s^ , in a white 
channel. 



1 Introduction 

In spite of lengthy and persistent investigations, at present we have no firm determination for 
the mass of the a meson - the resonance in the 280-900 MeV region. This resonance reveals 
itself in the vrvr channel as a pole in the complex-M plane, in the (JJ^*" = 00++) partial 
wave. Numerous calculations produced mass values distributed over all the low-energy interval 
^/s = M <90Q MeV, with various widths from 200 MeV up to 1000 MeV. Such a situation 
emerged in the nineties [1]. The results of the latest analyses are clustered in a smaller mass 
region 400-600 MeV: see, for example, [2] (552t?^6)-^(232lf^) MeV and [Sj (48^17) -i(255±10) 
MeV) and the review of Bugg [4] 472±30 - i (271 ±30) MeV. The solution of the Roy equation 
at low energies produced a smaller mass UltfMeV - i(272l?3) MeV 0. 
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We see three sources for emerging uncertainties in the analyses of the vrvr amphtude near 
the threshold: 

(i) a not sufficiently good determination of the 00+"'' amplitude above M = 900 MeV, 

(ii) uncertainties in the definition of the left-hand cut in the tttt amplitude and 

(iii) uncertainties in low-energy tttt — t- tttt data. 

In the present paper we analyzed in detail all these sources of uncertainties. The examples 
considered in the paper demonstrate that the results obtained for the low-energy amplitude 
depend strongly on the assumptions made in the analysis. 

The meson spectra in the (00"^"*") wave were fitted by our group using the i^'-matrix technique 
O [71 E]- This technique provides us with an opportunity to fit simultaneously several reactions 
(such as TTTT, KK, rji], etc), taking into account correctly analytical properties and unitarity 
in all investigated amplitudes. This way we have determined the resonance structure of the 
scalar-isoscalar wave at 500 < ^/s < 1950 MeV; our results were summarized in [9]. 

However, in the i^-matrix amplitude the left-hand cut owing to crossing channels is deter- 
mined ambiguously (note that t and u channel meson exchanges depend on couplings and form 
factors, which are not well known). The impossibility to write down precisely the contributions 
of left-hand cuts leads to a freedom in the interpretation of the tttt — t- tttt amplitude in the 
y/s < 500 MeV region. In our i^-matrix analyses [01 El E] of the isoscalar scalar wave we 
modeled the contribution from the left-hand cut at s < by introducing several poles in this 
region with fitted parameters. Describing this partial amplitude in the region 280 < ^/s < 1900 
MeV, we usually did not observe a pole which could be interpreted as the a-meson. However, 
in some solutions (not the best ones) such a pole appeared. 

Having this background, we fitted in [lO] the amplitude OO"*""*" in the region 280 < y/s < 900 
MeV separately in the framework of the dispersion relation approach sewing the N/ Z)-solution 
with the i^-matrix one at 450 < ^/s < 900 MeV. As a result, the best fit, accounting for the 
left-hand cut contribution (it was a fitting function), contained the cr- meson pole at Mo- = 
(430 ± 150) - i(320 ± 130) MeV [lOj. 

One can think that the ambiguity problem may be solved with the help of the investigation 
of the vrvr scattering in all three {u,d,s) channels (see [11] and references therein). However, 
this procedure requires the analytical continuation of the pole terms into regions being rather 
far from the pole mass. This supposes the knowledge not only of both resonance form factors 
and the energy dependence of resonance widths. The high spin states lead to the divergence in 
crossing channels. It is only the summing over all sets of states that resolves these divergences 
resulting finally in the Regge behavior and therefore, requires model-dependent calculations. 

The /^-matrix analysis O [71 E], being performed at a distance from the left-hand cut, gives 
masses and full widths of resonances (i.e. the position of poles) as well as the residues of the 
poles, namely, couplings of resonances to different channels. These couplings are factorized; 
this is a criterium for dealing just with a particle, though unstable. Besides, the coupling 
interrelations allow one to define the quark content of a particle, provided this is a qq state. 
This way the states found in the i^-matrix analysis can be classified as qq nonets. The K- 
matrix analysis determines two nonets and one extra state in the 600-2000 MeV region. One 
of the possible classifications is given in [9]: 

[/o(980), /o(1300)]„=i, [/o(1500), /o(1750)]„=2 ; 

where n = 1,2 are the radial quantum numbers. Here the broad state /o(1200 — 1600) and 
the cr-meson are superfluous for the qq nonet classification. The position of resonances in the 
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Figure 1: Complex-M plane for the (JJ^*^ = 00"'""'") mesons [9]. The dashed hne encircles 
the part of the plane where the ii"-matrix analysis [7j reconstructs the analytical i^-matrix 
amplitude: in this area the poles corresponding to resonances /o(980), /o(1300), /o(1500), 
/o(1750) and the broad state /o(1200 — 1600) are located. Beyond this area, in the low-mass 
region, the pole of the light cr-meson is located (shown by the point the position of pole, M = 
(430— z320) MeV, corresponds to the result of N/D analysis ; the crossed bars stand for a-meson 
pole found in |I0]). In the high-mass region one has resonances /o(2030), /o(2100), /o(2340), 
see [H]. Solid lines stand for the cuts related to the thresholds tttt, tttttttt, KK, rjrj, rj-q'. 



IJ^'-' =00^^ wave is shown in Fig. [H 

In this classification the broad state /o(1200 — 1600) is a glueball descendant [121 IlSj- Due 
to another classification the broad /o(1200 — 1600) state belongs to the first nonet and the extra 
state is /o(1300). Both these states are flavor blind and one of them is superfluous for the qq 
systematics. The /o(1200 — 1600) state acquired a large width because of the accumulation 
of widths of neighboring states: in nuclear physics such a phenomenon had been studied in 
[m [151 [IB] , in meson physics in . 

The paper is organized as follows. In Section 2 we provide formulae used in the i^-matrix and 
D-matrix approaches. In Section 3 we discuss an idea of the confinement singularity 1/s^. Such 
a singularity in the t-channel (1/t^) corresponds to the linear rising potential which describes 
meson spectra in qq [19j, bb [20] and cc [21j channels and gives correct values for the partial 
widths of radiative and hadronic decays of confined qq states [22j|. Although this singularity is 
expected to be in the color octet state, it can have also a color singlet component and appear 
in the s-channel. The ii'-matrix and D-matrix analyses of the 00"'""'" wave in the energy interval 
280 < ^/s < 1900 MeV are presented in Section 4. In the Conclusion we summarize the results 
concentrating on the low-energy region. 

Some clarifying points are made in the Appendices. In Appendix A the dynamical mech- 
anism of the singularity 1/s^ in the qq — )■ qq amplitude is discussed. A simple description 
of the low energy vrvr scattering (270 < ^/s < 900 MeV) in terms of the dispersion relation 
approach which allows to incorporate easily the singularity 1/s^ into the analytic and unitary 
amplitudes is given in Appendix B. In Appendix C we present the unitary tttt scattering in the 
threshold regions taking into account the mass differences of the vr+vr" and vr^vr" systems which 
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are essential for the extraction of a^. 

2 The i^-Matrix and i^-Matrix Techniques 

Here we discuss the analytic properties of amplitudes restored in terms of the i^-matrix and 
D-matrix techniques. 

2.1 The iC-matrix approach 

For the 5*- wave interaction in the isoscalar sector we use, as previously [7], the 5-channel K- 
matrix: 

where K^j^ is a 5x5 matrix {a,b = 1,2,3,4,5), with the following notations for meson states: 
1 = TTTT, 2 = KK, 3 = rjf], 4 = 7777' and 5 = multimeson states (four-pion state mainly at 
^/s < 1.6 GeV). The are coupling constants of the bare state a to meson channels; the 
parameters fab and Sq describe the smooth part of the i^'-matrix elements (sq > 1.5 GeV^). 
The factor (s — sa)/{s + sao), where sa ~ (0.1 — 0.5)m^, describes Adler's zero in the two-pion 
channel. However, in the K-matrix analysis we introduced this factor also in other channels to 
suppress the effect of the left-hand side false kinematic singularities in the i^-matrix amplitude. 



Figure 2: Graphical representation of the spectral integral equation for the K-matrix amplitude. 



2.1.1 Spectral integral equation for the K-matrix amplitude 

Discussing meson-meson scattering and production amplitudes, we use the dispersion relation 
(or spectral integral) technique. In terms of this technique we write for the K-matrix amplitude 
a spectral integral equation which is an analog of the Bethe-Salpeter equation [23] for the 
Feynman technique. The spectral integral equation for the transition amplitude from the 
channel a to channel b is presented graphically in Fig. |2] and reads: 

= / ^^^:^^PAs')KAs\s) + Kab{s) . (2) 

Here Pj{s') is the diagonal matrix of the phase volumes, Aaj{s, s') is the off-shell amplitude and 
Kjb{s, s') is the off-shell elementary interaction. Let us remind that in the dispersion relation 
technique, just as in quantum mechanics, there is no energy conservation for the intermediate 
states. For more details see [9], Chapter 3. 
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The standard way of the transformation of eq. ([2]) into the K-matrix form is the extraction 
of the imaginary and principal parts of the integral. The principal part has no singularities 
in the physical region and can be omitted (or taken into account by a re-normalization of the 
K-matrix parameters): 



/ 
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iAaj{s,s)pj{s)Kjb{s) (3) 



For the amplitude Aab{s) one obtains the standard K-matrix expression which in the matrix 
form reads: 

A = Aipk + k, or A = k{i-ipk)-\ (4) 

The factor [I — ipK)~^ describes the rescattering of mesons, it is inherent not only in two-meson 
transition amplitudes but in production amplitudes as well. The P-vector method describes the 
production of particles in cases when an initial interaction should be taken into account only 
once, for example, for the production of mesons from the 77 collision or from proton-antiproton 
annihilation: 

A,{pp) = pA{i - ipk)-^] (5) 
Elements of the vector Pj have a form similar to the K-matrix elements, eq. ([1]): 

A 0^"' 

The first term in eq.(l6]) refers to the production of resonances; the second one, Fj, to a non- 
resonant production. 

The standard form of the two-particle phase volume is 



Pa{s,mia,m2a) = ]l ^ , a = 1, 2, 3, 4 (7) 

where rriia and are masses of the final particles. In the case of different masses this 



expression includes the term — {rriia — rn2aY which in the K-matrix approach can be a 
source of false kinematic singularities on the first (physical) sheet: the loop diagram amplitude, 
B{s), does not contain this type of singularities. Such a cancelation can be taken into account 
effectively by replacing the rjrj' phase volume: 



{s - {mia + m2aY){s - {mia - m2aY) ^ s - {mia + m2a) 



S2 



For the restoration of the amplitude we need to take into account not only the cuts related to 
threshold singularities of the stable particles but non-stable ones as well. In the 00"'""'"-amplitude 
the four-pion state gives cuts related to 71717171 (at the real s-axis, ^/s = 4:m„) and in the complex- 
s plane related to the production of vector and scalar particles: nnp (at ^/s = 2m^ -|- mp with 
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a complex mass nip), pp (at = 2mp) and /o/o- Let us write the phase space factor for the 
pp-state which contains An, nnp and pp threshold singularities: 



dSi2 f dss4^ 2 



P47r(s) = / / Gi^{s,Su,Ss4) p{s,^/s^,^/s^) X 

J IT J IT 

4 m| 4 

G'^(si2)(si2-4m2)p(si2,m^,m^) (^^(534) (S34 - Aml)p{ssi,m^,m^) 



{si2 - Mjf + (M,r,)2 (S34 - Mjf + (M,r, 



(9) 



The form factors ^171(5, S12, S34), G'(si2)(si2), ^(534) are introduced into ^ to provide the 
convergency of the integrals. This phase volume describes production of pp in the S'-wave and 
P-wave production of pions in the p-meson decays. Being near a pole, hadronic production 
cuts split this pole into several ones located on different sheets of the complex-s plane. 



2.2 The D-matrix approach 

The considered above approaches allow us to distinguish between "bare" and "dressed" parti- 
cles: due to meson rescattering the bare particles, with poles on the real-s axis, are transformed 
into particles dressed by "coats" of meson states. In the K-matrix approach we deal with a 
"coat" formed by real particles - the contribution of virtual ones is included in the principal 
part of the loop diagram, -B(s), and is taken into account effectively by the re-normalization of 
mass and couplings. 

In the dispersion relation D-matrix approach one can take into account the "coat" of virtual 
mesons. The D-matrix amplitudes describe transitions of bare states. 

Let us consider the block Dap which describes a transition between the bare state a (but 
without the propagator of this state) and the bare state [3 (with the propagator of this state 
included). For such a block one can write the following equation: 

Dap = ^ Bi^ndriP + daj3 (10) 

Or, in the matrix form: 

t) = DBd + d D = d{I - Bd)-^ (11) 
Here the c? is a diagonal matrix of the propagators: 

d = d^ag . . . , P„ i?2 • • •) (12) 

where Ra are propagators for non- resonant transitions (discussed below), and the elements of 
the i?-matrix are equal to: 

R 'Pjjs ,mi,,m2j)gj' 

= l^B^, = 737—^ • (13) 

The gf''^ and gf^"'"' are right and left vertices for a transition from the bare state a to the 
channel j. For the pole terms there is a clear factorization: 
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However, non-resonant terms do not provide such a factorization. A solution of this problem 
is to introduce for non-resonant transitions a separate propagator and vertices from every initial 
state i. Moreover, for the description of the non- resonant terms between different initial and 
final states a second propagator with permutated left and right vertices is needed. In this 
case the propagator index provides automatically a unique identification of the transition term. 
Then for non-resonant transitions from the vrvr channel we have: 

9^^""^'^ Rigf^""^'' + g^^'^^'^R^gf''^'^ (15) 

where N is the number of pole terms. The non-zero left and right vertices can be taken as: 

■' S + Sq 

L(N+2) _ ^ RiN+2) _ r 1 GeV^ + So 

■' S + Sq 

and 

R(N+1) _ L(N+2) _ BiN+2) _ „ 

gj>i ~ gj>i — fi — u i-L'j 

Another alternative parametrization for the non-zero terms is: 

L(N+i) _ r R(N+i) _ . p _ 1 GeV^ + So 

gj - hj gi - J- Ri : , 

L(N+2)_^ R{N+2) _ r p _ 1 GcV^ + So , . 

gi - 1 gj>i - hj R2 : • lioj 

s -\- So 

With such a definition the amplitude Aab is the convolution of the matrix Dq,^ with right 
and left coupling vectors, g^^'°'^ and gi^'^^: 

Aa, = Y.ga''"'^do.o.D^m^^\ (19) 

The P-vector amplitude has the form: 

^ = E P^'^^d^^D^pg^^^^ P = (Ai, A2, . . . , A„, . . .) (20) 

where couplings A^ and non-resonant terms Fj are the same as in eq.([6]). 

In the present fits we calculate the elements of the 5;?^^ using one subtraction taken at the 
channel threshold Mj = (my + m2j): 

is) - (M') + (s- M') 7ds' gf''''p,{s',rm„m2,)gf'^^^ 

ml 

In the case of the non-resonant terms parameterized in the form f ll8l) and the S-wave vertices 
parameterized as constants the expression for elements of the 13 matrix can be rewritten as: 



oo 



V 



ds' Pj{s',mia,m2a) I L{/3) 



TT {s' - s - iO){s' - 



Table 1: Coupling constants given by quark combinatorics for {qq)i=o meson and glueball decays 
into two pseudoscalar mesons in the leading terms of the l/Nf. expansion. The $ is the mixing 
angle for nn = [uu + dd) / \/2 and ss states: nn cos $ + ss sin $. The is the mixing angle for 
7] — 7]' mesons: rj = nn cos B — ss sin 6 and rj' = nn sin B + ss cos B with B ^ 37° 
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g cos$/\/2 
g cos$/\/2 
g{y/2sm.^ + a/a cos $) 
51(^2 sin $ + v^cos$)/\/8 
5f(cos^B cos$/V2 + v^ sin$ sin^ B) 
51 sin B cosB(cos$/V2 - sin$) 



gg state decay coupling 



G 

G 
VXG 
^/XG 

^(cos^B + yAsin^B) 
G(l - A) cos B sin B 



identity 
factor 

1/2 

1 

1 

1 
1/2 

1 



where the parameters b' depend on decay channels only. 

In the case of the D-matrix approach it is not needed to introduce the regularization of the 
riT]' phase volume and, therefore, we use the standard expression ([7]). It is also not necessary 
to introduce any regularization for the D-matrix elements at s = 0: this point is not singular 
in this approach. Thus, in the D-matrix fits, the term with the Adler zero was introduced in 
the TTTT channel only. Technically, it can be done either by the modification of vertices or by 
the modification of the vrvr phase volume: 



/ . s-sa s-Am^ 

Pi(s, m^, m^) = — ■ W (23) 

s + Sao V s 

For qq states one can relate the decay couplings g^^ in terms of the rules of quark combina- 
torics (see [9j, Chapter 2, and references therein). The couplings for channels a = vrvr, KK, rji], 
Tjrj', calculated in leading terms of the 1/Nc expansion, are presented in Table [TJ The couplings 
depend on the constant g which is universal for all nonet states, the mixing angle $ which 
determines the proportion of the nn = {uu + dd) / \/2 and ss components in the decaying qq 
state, and the ss production suppression parameter A ~ 0.5 — 0.7. Two scalar-isoscalar states 
of the same nonet are orthogonal if: 

$(^) _ = ±90". (24) 

The equality of the coupling constants g and the fulfilment of the mixing angle relation flM|) is 
a basis for the determination of mesons of a gg-nonet. 

The gluonic states are decaying in the channels a = vrvr, KK, rji], rji]' with the same couplings 
as the gg-state but at a fixed mixing angle $ — ^giuebaii which is determined by the value of 



A, namely: ^giuebaii = cos ^ ^2/(2 + A). The corresponding couplings are given in Table [T]as 
well. 
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3 Confinement Interaction in the qq Sector 



The description of mesons of the qq sector is a source of information about quark confinement 
interaction. These interactions contain t-channel singularities of scalar and vector type. The 
t-channel exchange interaction can be both in white and colour states, c = 1 + 8 though, of 
course, the colour-octet interaction plays a dominant role in meson formation. 

The observed linearity of the gg-meson trajectories in the (n, M^) planes [21], where n is 
the radial quantum number of the gg-meson with mass M, provides us the t-channel singularity 
Vconf ~ l/?"^ or, in coordinate representation, Vconf ~ In the coordinate representation the 
confinement interaction can be written in the following potential form [9l [19] : 

Vconf = (/ ® /) 65 r + (7m ® 7m) bv r , (25) 
bs ~ -by ~ 0.15 GeV-^ . 

The first term in 0251) refers to scalar interaction (/ /), the second one to vector (7^ 7^) - 
in the qq sector the scalar and vector forces are approximately equal. 



3.1 White remnants of the confinement singularities 

We have serious reasons to suspect that the confinement singularities (the t-channel singularities 
in the scalar and vector states) have a complicated structure. In the color space these are 
octet states but, may be, they contain also white components. The octet exchange interaction 
contains quark-antiquark and gluonic blocks. Therefore, the question is whether V^^^j^^^^^^^lq'^) 

has the same singular behavior as vj:^l^j^^^^^^^{q^). The observed linearity of the (ra, M^)- 
trajectories, up to the large-mass region, M ~ 2000 — 2500 MeV [23], favors the idea of the 
universality in the behavior of potentials V;I,^i/i„emeni and V^^l^^^^^^^^ at large r, or small q. To 
see that, let us consider, as an example, the process 7* — )■ qq, Fig. |3^. We discuss the color 
neutralization mechanism of outgoing quarks as a breaking of the gluonic string by newly born 
gg-pairs, see the discussion in [2S]- At large distances, which correspond to the formation of 
states with large masses, several new gg-pairs should be formed. It is natural to suggest that a 
convolution of the quark-gluon combs governs the interaction forces of quarks at large distances, 
see Fig. [3)d. The mechanism of the formation of new gg-pairs to neutralize color charges does 
not have a selected color component. In this case all color components 3 ® 3 = 1 + 8 behave 
similarly, that is, at small g^ the singlet and octet components of the potential are uniformly 

singular, Ko2/memeni(?^) ~ VcolfinemenM^) ~ V?^- 

If the confinement singularities have, indeed, white constituents, this raises immediately the 
following questions: 

(i) How do these constituents reveal themselves in white channels? 

(ii) Can they be identified? 

In the scalar channel we face the problem of the a meson [IJ^'-' = 00^^): what is the nature of 
this state? If the white scalar confinement singularity exists, it would be reasonable to consider 
it as the a meson revealing itself: because of the transitions into the tctt state, the confinement 
singularity could move to the second sheet. If so, the a meson can certainly not reveal itself as 
a lonely amplitude singularity 1/t^ but a standard amplitude pole or a group of poles. 

A similar scenario may be valid also for the vector confinement singularity in the tttttt 
{IJ^'-^ = 01 ) channel. In this case it is natural to assume that the white confinement 
singularity couples with the channel pn, splits and dives into the complex- M^^^ plane. 
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Figure 3: a) Quark-gluonic comb produced by breaking a string by quarks flowing out in the 
process e~^e~ — )■ 7* — )• gg — )■ mesons, b) Convolution of the quark-gluonic combs, c) Example 
of diagrams describing interaction forces in the qq systems. 

An illustrative example of a set of loop diagrams of the Fig. [3t type is considered in 
Appendix A. In this example we demonstrate how the strong singularity, may arise in 

scalar and vector channels of the interaction block. An analysis of the simple case when we 
have only two poles on the lower part of the second s-plane is performed in Appendix B: the 
obtained result is in agreement with those obtained in analyses performed over low-energy data 
only t26l[271[28l[29]. 

4 The iT-Matrix and D-Matrix Approaches in Fits to 
the Data at 0.28 < < 1.95 GeV 

Here we present a comparative analysis of the results obtained with the i^'-matrix and D-matrix 
methods. These approaches give rather similar results for the /o-resonances at ^/s < 2 GeV. 
In Table [2] we show the data used in these analyses and give corresponding for different fits. 
In Table [3] we list the masses of bare states, mixing angles and other parameters used in the 
minimization procedure. 

4.1 The matrix fit 

In the analysis of the present data set we fitted data in two steps. In the first step all couplings 
were optimized as free parameters; in the second step we imposed relations Table [1] for the 
poles with masses above 1 GeV. We did not observe any deterioration of the data description 
due to these restrictions but a rather notable improvement in the convergency of the fits. For 
the lowest K-matrix pole we do not impose any constraints: the global coupling and mixing 
angle for this pole given in Table [3] are simply calculated from the couplings into the tttt and 
KK channels. 

In the present solutions there are two candidates for a glueball: it is either the third or the 
fourth K-matrix pole (with a mass around 1200 MeV). For the glueball candidate we introduced 
in addition a glueball decay coupling (see Table H]). However, this coupling provided only a 
small improvement and did not allow us to distinguish between these two cases. 

The fit is hardly sensitive to the tttttttt couplings for the two lowest K-matrix poles; in the 
final solution we fix them to be zero. 
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To get a combined description of all reactions, we introduced non-resonant terms for the 
transition from the tttt channel to other final states. We did not find a notable sensitivity to 
non-resonant transitions between other channels. 

The data can be described with a very small re-optimization of the K-matrix parameters 
found in [HI [TJ |8]. We did not find any change in the pole structure of the scalar-isoscalar 
amplitude above 900 MeV. However, one of the pole singularities situated around s = moved 
to higher masses. Its position, as well as the positions of other poles, is given below in f l26p . 

4.2 The D-matrix fits 

D-matrix parameters can be expressed in the same terms (bare masses and couplings) as pa- 
rameters of a K-matrix fit. The subtraction point for calculation of the real part of the loop 
diagrams, i.e. i?^^(M^) in eq. (12T|) . was taken at the corresponding two-particle threshold; the 
parameters were optimized in the fit. In such an approach our data base can be described 
with a very similar quality as in the framework of the K-matrix approach, see Table HI As ex- 
pected, the D-matrix fit provides a better description of the K^i data due to the more correct 
behaviour of the amplitude near left-hand side singularities. The behaviour of the phase shift 
6q and its description in the mass region from the threshold to 1 GeV is shown in Fig HJ 

Below we present four D-matrix solutions: the bare masses and their couplings are given 
in Table [3] (Solutions 2,3,4,5). In Solutions 2,3 the K^^ point near 500 MeV was taken with 
the error given by the experimental group. However, these solutions do not reproduce this 
point satisfactory. To force the tttt phase shift to describe this point, we decreased the error 
by a factor 10 and repeated the D-matrix fit of the data. In such an approach we were able to 
describe the data at 500 MeV rather well (Solutions 4,5); however, we obtained a systematically 
worse description of the proton-antiproton annihilation into the vr^vr^vr'^ and r^r^vr^ channels (see 
Table [2]). 

The solution with the 1/s^ term included (Solutions 3,5) produced a better total ^ 
slightly better description of the Kg^ data. The term 1/s^ can produce two additional poles 
in the mass region below the tttt threshold. The pole in the mass region around 400 MeV has 
moved to lower masses by about 80 MeV compared to solutions without the 1/s^ term, see 
(|26|) . while the poles situated above 900 MeV practically do not change their positions. 

It is seen from Table [3] that the masses of bare states are hardly changed from the K-matrix 
solution and most of the couplings are shifted by less than 20%. The positions of the amplitude 
poles above 900 MeV also changed very little: 





Solution 1 


Solution 2 


Solution 3 


Solution 4 


Solution 5 


(T-meson 
/o(980) 
/o(1300) 
/o(1500) 
/o(1750) 


420-i 395 
1014-i 31 
1302-1180 
1487-i 58 
1738-i 152 


407-1281 
1015-i 36 
1307-1137 
1487-i 60 
1781-1140 


365-i 283 
1012-i 31 
1303-i 140 
1483-i 55 
1787-i 143 


414-i 186 
1005-i 20 
1332-i 140 
1487-i 55 
1795-i 109 


406-i 192 
1005-i 23 
1326-1137 
1486-i 55 
1794-i 114 



(26) 



The relative position of the poles and the threshold singularity cuts is demonstrated in Fig. [TJ 

We see that a fit of the Kf>4, data with the use of the D-matrix approach unambiguously 
reveals the pole in the mass region around 300-400 MeV, the low-mass a-meson. 

To trace the origin of the a pole, we multiplied all couplings by the factor /3 and the 
non-resonant terms by and scanned this parameter from 1 (the physical amplitude) to 
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<s, MeV <s, MeV 



Figure 4: Description of the i^'e4 data with the /^-matrix solutions 2 ,3 (with standard errors 
for the point 6Q{500MeV) |30]) and solutions 4,5 (with the decreased error for this point). 



(amplitude with poles corresponding to the bare masses). Such an investigation shows that the 
(j-pole is originated from the Adler regularization term. In the best fit the Adler regularization 
point is optimized rather close to the physical region sao — 0.15 GeV^. To check the stability 
of this point we have performed the fit with this point fixed at sao = 0.5 sao = 1 and s^o = 1-5 
GeV^. We observe a small deterioration of the total due to a worse description of the 
7171 —7- KK and tttt — )■ rjr] amplitudes. However, the fit with sao = 0.5 GeV^ gives the best 
description for the proton- antiproton annihilation into the tt'^tt^tt^ channel that is one of the 
most sensitive reactions to the description of the low tttt mass region. The positions of the 
poles in all three solutions coincide remarkably well and hence, we conclude that the position 
of the cr-meson depends very little on the exact position of the Adler regularization term. 

4.3 Calculation of the scattering length 

In our expression for the tttt scattering amplitude which takes into account the and TT^TT 
phase volumes the vrvr-phase does not go to zero on the threshold of two charged pions, see 
Appendix C. We calculate the scattering length of the tt+tt" system at the threshold of two 
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charged pions using the following expression: 



,1/2 



-Re 



• s:(0) iS^°> 

sm6n 'e 



with ao = ^a^(t\ 



fe-s>0 
1 



k = -Js 



(m^++m^-)2 and m,, 



The scattering length values extracted from the D-matrix solutions are equal to: 



Solution 2 


Solution 3 


Solution 4 


Solution 5 


0.253m;i 


0.209m;i 


0.204m-i 


0.177m;i 



(27) 



(2^ 



It is seen that the inclusion of the term decreases the scattering length by ~ 0.05 m^^ and 
a similar effect comes from a precise description of the i^'e4 experimental point at 500 MeV. 

The amplitude phase was extracted by the K^^ collaboration under the assumption that 
it is equal to zero at the threshold of two charged pions. Then there is a question about 
the uncertainty which appears when these data are fitted with an expression which takes into 
account exactly the thresholds of neutral and charged pions. To check this we put in Solution 
2 all pion masses equal to the mass of a charged pion. As expected, notable deteriorations 
were observed only in the proton-antiproton annihilation into three neutral pions and at low 
energy points for the Kf,^ data. With a very small tuning of the parameters we obtained very 
similar values for the description of the K^^ data. The scattering length which in this case 



is calculated as ao 



3.(±) 
2"'0 



appeared to be 0.248 m^^. Then, with this parameters fixed, we 



introduced back the difference between neutral and charged pion thresholds but not refitted 

a value 



-1. 



the data. The scattering length obtained by eq. ( !27l) was found to be 0.260 ml 
which is very close to that obtained in Solution 2. Thus, we conclude that the investigated 
uncertainty is less than 0.010 and is smaller than the systematic error which comes from 
different parameterizations of the amplitude. 

It is instructive to compare the results of eq. (128|) with those obtained without taking into 
account different values of 7r°7r° and 7i~^7i~ threshold singularities: ao = (0.233 ± 0.013)/x~^[3j, 
ao = (0.220 ±0.005)/i;^^| 



5 Conclusion 

The analysis of the large data sets performed in the framework of the i^-matrix and D-matrix 
approaches demonstrates a very good stability for the amplitude parameters and pole positions 
above 900 MeV. The description of the Ke4 data demands the presence of the pole slightly 
above the tttt threshold. The pole position was found to be at 390 ± 45 — z 210 ± 50 MeV and 
the scattering length 0.215 ± 0.040 m~^. 

The confinement singularity, 1/s^, slightly improves the overall description but is not crucial 
for a good description of the Kg4^ data and for the existence of the a meson pole singularity. 
However, the presence of such a term influences the pole position shifting it to lower masses by 
about 100 MeV and shifting the scattering length to lower values (by ~ 0.05 m'^). 

The imaginary part of the a pole position in the solutions which fit precisely the data point 
at 500 MeV is lower by about 100 MeV compare to the solutions where the fit is not forced 
to describe this point. The scattering length in such solutions is also systematically shifted to 
lower values by ~ 0.05 /i^^ 
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Within the description of the OO+^-wave in the channels vrvr, vrvrvrvr, KK , rjr]' we obtain the 
following complex masses of the /o resonances: 

/o(980) M/ = 1011±5 -i31±4 MeV 

Mil = 1035±50 - i460ib50 MeV 
/o(1300) M = 1285±30-il60±20 MeV 

/o(1500) M = 1488±4 -i53±5 MeV 

/o(1790) M = 1775±25-il40±15 MeV (29) 

The masses of the D-matrix approach, eq. fl29p . coincide well with those obtained in the K- 
matrix approximation [9j. The /o(980) is determined by two poles, on the second (under the 
TTTT threshold) and third (under the vrvr and KK thresholds) sheets - the same splitting of poles 
we have in the i^'-matrix solutions |32j . 

For the low mass region the solution with singularity gives several poles on the second 
sheet: 

fo{ai) M = 365±15-i283±12MeV 
fo{an) M = 80±10-il87±15MeV 
foiaiii) M = -94±12-i93±10 MeV. (30) 

If the fit is forced to describe the Kf.4 experimental point at 500 MeV, we have: 

/o(a/) M = 406±15-U92±15 MeV 
/o(cT//) M = 74±10-il90±50 MeV 
foiajjj) M = -96±22 -il00±25 MeV. (31) 

We also test the changes in the description of data with an elimination of the 1/s^ singularity. 
In this case the fit to the data gives the masses of the /o resonances at y/s > 900 MeV practically 
the same as in ref. (129|) - the changes are in the low-mass pole structure. Without the 1/s^ 
singularity, the position of the a pole in the fit, neglecting the 500 MeV point, gives: 

foiai) M = 407±12-i289±10 MeV 

(32) 

and with the fit forced to describe the 500 MeV point: 

/o(a/) M = 412±12-U86±15 MeV (33) 

So, the (j-meson arises as a pole near the vrvr threshold in both versions, with and without 
including the confinement singularity (1/s^) into the tttt scattering block. Though the confine- 
ment singularity leads to the appearance of several poles under the tttt cut, it is hardly possible 
to distinguish these two versions on the basis of the data. 
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a) 




b) 




Figure 5: Interaction block (a) and sets of loop diagrams for S and V exchanges. 



6 Appendix A: Examples of Sets of Diagrams Resulting 
in Singularities 

Here we consider, as an example, the confinement set of the loop diagrams. Fig. [3b, and present 
an illustrative calculation which results in singularities of the type in scalar and vector 
channels. We use the following interaction blocks, see Fig. [5^.: 



+G^v^ ( {h - hf U{h)^,i^{k,)G'^' ( {k[ - k',f UiKhMK) ■ 



(R) 



(34) 



Then the confinement interaction turns into a set of the loop diagrams. Fig. [5h> and Fig. 
The scalar and vector exchanges, correspondingly, read: 



y(,. Bsjt) By it) 

For a scalar loop diagram one has: 

4m2 



d^2iP';k[,-k',: 



^' 2(27r)2 
Sp[(k'2 - m)CK + m)] ^2{t' - Am^ 



^ d%d%S{k'^ - m^)6{k'^ - m^) ^ ^ 



16n 



Am?' 



(35) 



(36) 



Here we replace G^^^G^^^ — )■ Ng. An analogous loop for vector exchange reads: 
-5^^By{t)= J —d<l>2{P;k^,-k2) 1' -t- iO 

4m2 



^ fit/ 



dt' Gl{t'){2m^ + t') 1 



Am? 



TX 



Am? 



t' -t-iO 



IGtt 
15 



(37) 



Confinement singularities appear if Bs{t) and By{t) behave near t = as follows: 



Bsit) = 1 - 
that means the requirements 

dt 



(3s 



Bsit) 



Oif), Byit) = 1 



/3v 



i=0 



dt 



Bv{t) 



0. 



t=o 



(3^ 



(39) 



7 Appendix B: Simplified Consideration of 
the 00++ Wave in the Low-Energy Region 

The partial pion-pion scattering amplitude being a function of the invariant energy squared, 
s = M^, can be represented as a ratio N{s)/D{s) jBSj), where N{s) has a left-hand cut due 
to the "forces" (the interactions caused by the t- and w-channel exchanges), and the function 
D{s) is determined by the rescattering in the s-channel. The standard presentation of the 
N/D- method may be found, for example, in pM]). 

The 7171 scattering block related to the 1/s^ singularity reads: 



G{s)-G{s) 



(40) 



The s-channel re-scatterings give a set of divergent terms which convolutes into the following 
unitary amplitude: 



A{s) = G{s)\g{s) + G{s)\u{s)\g{s) + ... 



G^(s) 



G'(s) 



" ds' G^{s')p{s') 



7T 



(41) 



4mI 



Here p(s) is the invariant vrvr phase space. In the physical region, at s > 4m^ and s on the 
upper edge of the threshold cut, we have: 



His) 



ds' G\s')p{s') 



TT S — S 



iO 



P 



ds' G^{s')p{s') 



71 



+ tG\s)p{s) 



(42) 



with the following relation to the JJ^'" = 00++ phase shift: p{s)A{s) = exp ^Z(5o(s)j sin5Q(s) 

The product of the vertices G^(s) is actually an A^- function, and we re- write G^(s) — )■ N{s)] 
this allows to present the amplitude fHTj) as 



Ms) 



N{s) 
D{s) 



D(s) 



4m,?. 



ds' p{s')N{s') 
n s' — s — iO 



The A'"-function, being determined by the left-hand singularities caused by forces due to t- 
channel and w-channel meson exchanges, is written as an integral along the left cut as follows: 



N{s) 



ds' L(s') 



7C S' 



(43) 
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where the value sl marks the beginning of the left-hand cut. For example, for the one- meson 
exchange diagram g'^/{m? — t) the left-hand cut starts at sl = 4m^ — m^, and the A^-function 
in this point has a logarithmic singularity; for the two-pion exchange, sl = 0. 

We replace the left-hand integral for A^(s), eg. (1431) . by the following sum: 

N{s) = J ^ 16. E , (44) 

— oo 

where L„ and s„ are "force parameters", — oo < s„ < sl- 

The pole approximation ansatz f l44p allows us calculate the scattering amplitude in the 
physical region: 



exp ( ^5',{s) ) sin 5',{s) = ' " . . (45) 

■s^ - E \/4m2 - Sn + iJs - 4m2 L„(s - s„) ^ 



Here we give an example of a very simple, and formally correct, consideration of the 00"'""'" 
wave in the low-energy region. Using eq. fH51) we write an analytic and unitary amplitude as 
follows: 

AUs) = J ^- . , (46) 

s2 - U/ + his + ig'^^s - 4/i2 j 

and hence, 

. 0/ _ -P^(g) _ [fc - (g + ib)] [k - {-a + ib)] [k-{c + id)] [k - (-c + id)] 
exp[ tdo{s}] - ^^^^^^ - ^^^^^^ _ _ ^^^^^^ ^^^^^^ _ _ ^^^^ 

with a > 0, b>0. (47) 

Considering (a, 6, c, 6) as parameters, we fit the data for 5q(s) in the energy interval 280 < 
^/s < 950 MeV, see Fig. El We obtain the following parameters and amplitude pole positions. 
Mi and Mif. 



Fig.4a : 


a = 


3.1/i^, 


b = LOyU^r, c 


= 7.7/i^, d = 


9.0/i^, 




Mi = 


(896.3 - 


i274.3)MeV, 


Mil = (2163.7 


-i2511.0)MeV 


Fig.4b : 


a = 


2.8/i,, 


,6 = 1.5/i^, , 


c = 7.7/i7r, d 


= 4.5/x^, 




Mi = 


(828.8 - 


i391.4)MeV, 


Mil = (2166.6 


-il252.1)MeV 


Fig.4c : 


a = 


2.6/i^, 


b = 1.9/i7r, c 


= hAjji^^, d = 


l.OyU^, 




Mi = 


(759.7- 


i509.1)MeV, 


Mil = (1529.8 


- 2275.4)MeV 


Fig.4d : 


a = 


2.5/i^, 


b = 1.3/x^, c 


= 50.0/i7r, (i = 


= 2.0/i^, 




Mi = 


(742.1 - 


i352.2)MeV, 


Mil = (14002.1 


3 - z559.9)MeV 



(48) 

In all solutions the scalar-isoscalar scattering length is not small: ~ (0.3 — 0.4)/i~^ . 

We use in the fit of Fig. |6] the values for 6q found in |6j in order to perform a more visual 
comparison of the obtained here results, eq. ( HHj) . with those in '^U\. Let us recall that we fit in 
[To] the amplitude 00^"'" in the region 280 < a/s < 900 MeV in the framework of the dispersion 
relation approach sewing the A^/D-solution with the K-matrix one at 450 < i/i < 1950 MeV. 
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Figure 6: Examples of the fit of low energy data [6l [35j in terms of eq. fj47|) 

Taking into account the left-hand cut contribution (it was a fitting function), we obtained in 
[ID] the best fit with the cx-meson pole at = (430 ± 150) - i(320 ± 130) MeV. So, the 
accounting for the left-hand cut and data at ^/s > 900 MeV results in a smaller value of the 

In the approaches, which take into account the left-hand cut as a contribution of some 
known meson exchanges, the pole positions were obtained at low masses as well. For example, 
the dispersion relation approach results: ~ (470 — z460) MeV [36], M^- — (450 — z375) MeV 
[37] , and the meson exchange models give: Mo- ~ (460 — iA50) MeV [38] , ~ (400 — iQO) 
MeV [39]. 



8 Appendix C: The nn Scattering Amplitude near Two- 
Pion Thresholds 

Here we consider the tttt scattering amplitude near two-pion thresholds taking into account the 
mass difference of charged and neutral pion systems, tt+tt^ and tt^tt^. 

The following vrvr-amplitudes describe scattering reactions near the thresholds: 



71 71 — )• 71 71 : A_ 



+ zA:0[(al° 
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1 xk_(x "^/lqCIqq ~1~ k^k_ AqqCI^^ ~1~ ('^— o)^] 



With kt = ^--fil,^K A;° = -^--/i^o = -v/FTA^. (49) 

Here = — /i^o — 0.07/i^+. The factor 1/2 in k^ arises due to the identity of pions in the 
7r°7r° state. 

We impose on the scattering length values the standard isotopic relations: 

2 1 

att = -ao{s) + -02(5), 

2 2 
«-o = -3«o(s) + -02(5), 



a, 



00 _ o^++ 

00 



2 _ 4 



2al^ + ai° = -ao(s) + -a2(s). (50) 



3 ^ ' 3 



Then at large k"^, when k"^ » A^, the unitary amplitudes of eq. fj49l) obey the isotopic relations: 

|ao(s) ^02(5) 



Alt 



1 — zfcao(s) 1 — 2/^02(5) 
|ao(s) §02(5) 



~° 1 — ikaQ{s) 1 — ika2{s)' 

aoo ^ l^o(^) I l^2(g) , . 

°° l-zfcao(s) l-^fca2(s)' ^ ^ 

The (/ = 0)-amplitude and the corresponding S'-matrix read: 

____ _ 2A__ - -A,, - A__ - -A_o, 

A++_lA+o 4++ _ 14+0 94++_l/ioo 

exp[ZiOo[s)\ - ^^^^ _ ^^^0^^ - ^2^++ _ ^^^0)* - (^2Att - |AgO)* • ' 

In the i^"^ — )■ e+z/(7r+7r^) decay the S'-wave pions are / = states, and the amplitude can be 
written as follows: 



Ai^K+ ^ e+z/(7r+7r-),=o,5-«,a.ej = A[l - tk',A',t + ik+Att] = 
1 _ ,-juo 00 _ jUO 0+ 

^1 - tklatt - iky^^o + k^okt[-a'o>-- + ("D'] 

Here the first term, A, is a direct production amplitude while the second and third terms take 
into account pion rescatterings. 

At large pion relative momentum, when fc^ >> A^, we have: 

a(k+ e+u{n+n-)j=o,s-wave) = A- ^— - . (54) 

V /fc2»A2 l-tkao{s) 

Recall that the factor (1 — ikao^s))"^ is due to rescatterings of pions in the 1 = state. 
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Tabic 2: 

List of the reactions and values for the i^-matrix and /^-matrix solutions: Solutions 3,5 



with taken into account confinement interaction, solutions 1,2,4 without it. 





Sol. 1 


Sol. 2 


Sol. 3 


Sol. 4 


Sol. 5 


N of 




K-matrix 


D-matrix 


D-matrix 


D-matrix 


D-matrix 


points 










~ 1/S2 





- 1/S2 






Tl 


le Crystal Barrel data 


from liquid H2: 














pp 7r°7r°7r° 


1.32 


1.37 


1.39 


1.44 


1.45 


7110 


pp n^rir] 


1.33 


1.34 


1.34 


1.33 


1.33 


3595 


pp 71^71^1] 


1.24 


1.33 


1.33 


1.55 


1.55 


3475 


from gaseous H2: 














pp 7r°7r°7r° 


1.39 


1.44 


1.45 


1.48 


1.49 


4891 


pp — >• n^rir] 


1.31 


1.34 


1.30 


1.43 


1.31 


1182 


pp 7r°7r°?7 


1.20 


1.22 


1.22 


1.31 


1.32 


3631 


from liquid H2: 














pp TT+TT^TT" 


1.54 


1.46 


1.45 


1.46 


1.47 


1334 


from liquid D2: 














pn — >■ 7r°7r°7r~ 


1.51 


1.47 


1.47 


1.46 


1.46 


825 


pn 7l~7l~7l~^ 


1.61 


1.54 


1.55 


1.50 


1.51 


823 


from liquid 7^2- 














pp KsKsTT^ 


1.09 


1.10 


1.10 


1.10 


1.10 


394 


pn K+K-TT^ 


0.98 


1.00 


1.00 


1.03 


1.02 


521 


pn KlK^tt^ 


0.78 


0.79 


0.79 


0.79 


0.79 


737 


from liquid D2: 














pp KsKsTT' 


1.66 


1.64 


1.64 


1.64 


1.63 


396 


pn KsK~Ti^ 


1.33 


1.31 


1.31 


1.31 


1.31 


378 




The GAM 


S data 




1.23 


1.13 


1.15 


1.32 


1.30 


68 


nn {r]ri)s-wave 


1.02 


1.05 


1.03 


1.58 


1.43 


15 


TTTT {m')s-wave 


0.45 


0.30 


0.35 


0.35 


0.34 


9 




The BNL data 


TTTT {KK)s-wave 


1.32 


1.13 


1.14 


0.97 


1.07 


35 






The CERN-Munich data: . 






TV 7r~^ — y TT TT 


1.82 


1.86 


1.86 


2.05 


2.03 


705 




The Ke4 decay data 


(5g(7r~7r+ 7r~7r+) 


1.51 


1.02 


0.84 


0.80 


0.83 


17 
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Table 3: 

The /Q'^''*^-resonances: masses Mn (in MeV units), decay coupling constants Qn of Table 1 (in 
GeV units), mixing angles (in degrees), background terms /„ and confinement singularity 

term G/s^ (factor G in GeV units). 





Solution 1 


Solution 2 


Solution 3 


Solution 4 


Solution 5 


Ml 


671 


685 


697 


611 


615 


M2 


1205 


1135 


1135 


1078 


1096 


M3 


1560 


1561 


1558 


1575 


1572 


M4 


1210 


1290 


1284 


1334 


1330 


Ms 


1816 


1850 


1848 


1858 


1857 


91 


0.860 


0.926 


0.892 


1.090 


1.083 


92 


0.956 


0.950 


0.935 


0.099 


1.066 


9s 


0.373 


0.290 


0.284 


0.302 


0.302 


9a 


0.447 


0.307 


0.308 


0.264 


0.275 


95 


0.458 


0.369 


0.370 


0.317 


0.330 


nil) 


-0.382 


-0.213 


-0.232 


-0.176 


-0.193 


^ 'in 


-0.322 


-0.500 


-0.500 


-0.500 


-0.500 


(1) (2) 

















(3) 

9J 


0.638 


0.534 


0.530 


0.511 


0.514 


(4) 

9471 


0.997 


0.790 


0.794 


0.691 


0.702 


(5) 
94J 


-0.901 


-0.862 


-0.856 


-0.797 


-0.814 


$1 


-74 


-83 


-82 


-81 


-82 


$2 


6 


-2.6 


-1.9 


-1.1 


-2.4 


$3 


9 


5 


5 


5 


5 


$4 


38 


31 


32 


25 


25 


$5 


-64 


-71 


-68 


-77 


-77 


/tttt— >7r7r 


0.337 


0.408 


0.358 


0.763 


0.687 




0.212 


0.036 


0.044 


0.103 


0.065 




-0.199 


-0.101 


-0.092 


-0.051 


-0.062 


fmr—^Tjri 


0.389 


0.438 


0.413 


0.538 


0.512 


f-K-K^rjrj' 


0.394 


0.518 


0.485 


0.610 


0.597 


G/s^ 








-0.00077/^2 





-0.00071/5^ 
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